The observed spectrum of a supernova remnant (SNR) is 
a superposition of many "local" spectra emitted by regions of 
SNRs that are under different physical conditions. The question 
remains as to whether the broadening of the high-energy end 
of the observed nonthermal spectrum of SNRs, like in G347.3- 
0.5 and SN 1006, can be an artifact of observations or it is 
a consequence of the microphysics involved in the accelera- 
tion process. In this note we study the influence of parame- 
ters variations (inside the volume and over the surface of SNR) 
on the shape of the high-energy end of the synchrotron (and 
also inverse Compton) spectrum. We consider three possibili- 
ties for these parameter variations: i) gradients downstream of 
the shock with constant maximum energy of the accelerated 
\Q electrons and the potential variation in time of the injection ef- 
CD ficiency, ii) then we add the possibility of the maximum en- 

■ ergy depending on time, and finally iii) the possible obliquity 
\ dependences of maximum energy and injection efficiency. It 

"q ■ is shown that gradients of density and magnetic field strength 
/"■n ' downstream of the shock are ineffective in modifying the shape 

■ of the synchrotron spectrum, even if an SNR evolves in the 
nonuniform interstellar medium and/or the injection efficiency 
varies in time. The time dependence of the maximum energy of 

| the electrons accelerated by the shock is also not able to make 

■ the observed spectrum much broader. The only possibility of 
r*"** | producing considerable broadening in the spectrum is the vari- 

■ ation in the maximum energy of electrons over the surface of 
SNR. In such a case, the obliquity dependence of the injection 
efficiency also affects the shape of the spectrum, but its role is 

\Q ' less significant. 
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Abstract. 



1. Introduction 

The properties of observed nonthermal emission from a num- 
ber of supernova remnants (SNRs) in X- and y-rays are the 
subject of many papers in the past decade because this emis- 
sion helps for understanding a number of important questions 
in modern high-energy physics and astrophysics (see Jones et 
al. rrWgl Drury et al. EOTfll Revnolds l2"U03l for reviews). 

Electrons accelerated by SNR shocks up to the maximum 
energies E max ~ 10 13 - 10 15 eV radiate their energy through the 
synchrotron emission in X-rays and through inverse Compton 
or nonthermal bremstrahlung in y-rays. The empirical expres- 
sion for the power-law spectrum of accelerated electrons with 
the upper cut-off (e.g., Gaisser, Protheroe, & Stanev 1998 1 



N(E) = KE S exp 



(1) 



is robust. The value a = 1 produces the most rapid cutoff in 
comparison with any physical model considered (Reynolds & 
Keohane 1999 and references therein) and, therefore, allows 
an estimate of the upper limit of E max from the X-ray spectra 
of SNRs (Reynolds & Keohane 1999 Hendrick & Reynolds 
I200TY 

The fit of the observed high-energy spectrum (X-ray and/or 
y-ray) of the given SNR may, however, require the broadening 
of the upper cut-off to the electron spectrum; i.e. the value of a 
would be less than unity and the high-energy end of the emis- 
sion spectrum of SN 1006 is broader than could be expected 
in the case a — 1 (Reynolds 1996 1. Detailed calculations give 
a * 0.5 for this SNR (Ellison et al.[2000l. The case of SNR 
G347.3-0.5 also requires the value a a 0.5 (Ellison et al. 1200 II 
Uchiyama et al. 120031 Lazendic et al. EOfffl . 

Such broadening can be caused by the microphysics in- 
volved in the electron acceleration process, such as in the 
model of Ellison et al 



BOOOt . On the other hand, this broad- 
ening could be an artifact of observations, where the observed 
spectrum is coming from a certain area of the projection of 



SNR on the plane of the sky. There is a part of the SNR surface, 
as well as a part of the SNR volume, projected onto this area. 
Different fluid elements are under different physical conditions 
in different places on the surface and in the volume involved. 
Therefore, the spectrum observed is a superposition of differ- 
ent "local" spectra. This "artificial" emission spectrum, which 
is a sum of the spectra produced under different conditions, 
is expected to differ from the spectrum radiated by a popu- 
lation of electrons characterized by a single set of magneto- 
hydrodynamic parameters. In this way, the "artificial" broad- 
ening of the electron distribution might be a consequence of 
the hydro-magnetic parameters variation in space, such as the 
suggested variations in E max (Reynolds 1996 1. Note, that we 
say "broadening" while not knowing in advance whether such 
an effect leads to broadening or narrowing of the spectrum end. 

The measure of broadening (a) in such a case depends on a 
number of factors. As an example, one of them could be a gra- 
dient of density downstream of the shock because the shapes 
of nonthermal spectra are sensitive to variations in the mag- 
netic field strength B and the normalization K inside the vol- 
ume; both of them are functions of density n whose profile, in 
turn, depends on the density gradient of the interstellar medium 
(ISM) before the shock. In this note we study the influence 
of the nonuniformity of the ISM, as well as other factors, on 
the artificial broadening of the high-energy end of the electron 
spectrum in the adiabatic SNRs. 

2. Models and results 

The spectral distribution of energy radiated by a single electron 
with energy E is 



p(E,v) = 



V3e 3 B D 



■F\- 



(2) 



where v is frequency, B p = -\/2/3B is an angle-averaged mag- 
netic field, F and v c oc BE 2 are the special function and the 
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Fig. 1. Influence of the broadening parameter a on the shape of 
the synchrotron spectrum (|6j. The values of a are marked near 
the respective lines, s = 2. 



critical frequency known from the synchrotron emission the- 
ory, e and m are electron charge and mass, and c is the velocity 
of light. The synchrotron emissivity of the population of elec- 
trons is 



P(v) 



= I N(E)p(E,v)dE 
Jo 



(3) 



and the synchrotron spectrum from SNR is given by the inte- 
gration over the given volume: 



S(v) 



-I 



P(v)dV. 



(4) 



Uniform plasma with the electron population having a 
broader distribution Q produces a "fit"-spectrum that does not 
require the volume integration 



Sfit(v) <* P(v). 



It can be written in the dimensionless form 

Xoo 



(5) 



(6) 



where e = E/E m . dx , e = v/v c (£' max , B), and the classic s 
2 is used. In these notations, most of the photons emitted by 
electrons with energy e have energy e m = 0.29e c , where the 
critical energy is e c = e 2 . Figure^demonstrates the role of the 
broadening parameter a in the modification of the shape the 
synchrotron spectrum. 

In the following sections, we calculate the high-energy end 
of the synchrotron spectrum from a part of SNR assuming that 
a = 1 and taking into account the possible variations of param- 
eters within the given volume (Sect. l2.lt or on the given surface 
(Sect. l2~3t of SNR in order to see whether we can obtain an ar- 
tificial spectrum with a shape similar to the fit-spectrum with 
a = 0.5. 



2.1. Variation of the parameters inside the SNR 
volume 

Let us consider a spherical SNR on the adiabatic phase of evo- 
lution in the ISM with the power-law density variation p (R) = 
AR~ W , where A and u> are constant; indexes "o" and "s" refer to 
the pre- and post-shock values. The dynamics of the adiabatic 
shock in this medium is given by the Sedov ( 1959 1 solutions 
where the shock velocity V oc /j-(3-w)/2 and R oc t 2 ^ 5 ^. These 
solutions are self-similar, which means that the distribution of 
parameters behind the shock, e.g. density n(a, t), can be written 
as n(a, t) = n s (f)«(«), where a is the Lagrangian coordinate, and 
the overline refers to normalized parameters, i.e. ones divided 
on their own values at the shock front, e.g. a = a/R. 

Let us assume, following most of the recent models of 
SNRs, that the ambient magnetic field is uniform. In the de- 
scription of the magnetic-field evolution downstream we fol- 
low the model of Reynolds & Chevalier ( 1981 1 and Reynolds 
J1998l >. where the strength of the magnetic field is B 2 = B 2 ± +B 2 
with the perpendicular and parallel components evolve inde- 
pendently behind the shock: B ± = nr/a, B\\ = (a/r) 2 . In this 
case, the behaviour of the magnetic field is also self-similar: 
B(a,& ) = B o cr B (0 o )B(a, O ) where the magnetic field com- 
pression factor ct b = ^1 + 16 tan 2 O ) / (l + tan 2 O )) 

We also follow Reynolds ( 1998 1 in the description of the 
evolution of electron distribution. Fluid element a = R(t[) was 
shocked at time t\. At that time the electron distribution on the 
shock was 



(7) 



(the value s — 2 is used throughout this paper). 

An electron loses its energy downstream due to the adia- 
batic expansion and radiative losses (synchrotron and inverse 
Compton on a cosmic microwave background). Its energy 
varies as (Reynolds 1998) 

E = £i«(a) 1/3 (l + I(a)Ei/E f y l (8) 
where I is the dimensionless function 

I(a, O , co, d) = ^ J' x^ /2 B eS (l) 2 n(l) in dx, (9) 

d = Bcmb/B , Bcmb = 3.27 pG is the magnetic field strength 
with energy density equal to that in CMB, and 



Et = 637/(B 2 ff s f)erg = \3B 5 2 t 4 l TeV 



(10) 



is the fiducial energy, B 2 ff = B 2 + 5c MB , B5 = 10 5 B e ff,s, U = 
10 

The conservation law for the number of particles per unit 
volume per unit energy interval 

a 2 dadE\ 



N(E,a,i) = N(Ei,a,tO- 



(11) 



4r 2 drdE' 

together with the continuity equation p {a)a 2 da = p(a, t)r 2 dr 
and the derivative dEJdE = h l ^E 2 /E 2 , implies that down- 
stream 



ME, a, t) = KE' 2 exp 



(12) 



4 
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with 

K(a,t) = K s (ti) n 4/3 a". 



(13) 



If K s oc V(t) b , the evolution of K is also self-similar down- 
stream: 



K(a) = K(a,t)/K s (t) = a [ib+uj{1 - b))l2 n(af\ 



(14) 



It is unknown how the injection efficiency (on which K de- 



pends) evolves in time. Reynolds ( 1998 1 considered three typi- 
cal alternatives for b as a free parameter, namely, b — 0, -1, -2 
(e.g. if acceleration efficiency is a constant fraction of post- 
shock pressure, then K oc V 2 ). Petruk & Bandiera (2006) show 
that one can expect b > and its value can even be as high as 
b * 5. 

The dependence E\(E) used in fl!2i can be obtained trivially 
from Eq. 0. It may be written in the form E-JE = S(a, E)~ l , 
where 



S = n(a) 1/3 -I(a,® )E/E f 



(15) 



In order to see the effect of the superposition of spectra 
from different regions inside the volume of SNR on the shape 
of the summary spectrum, we perform an integration along the 
radius of a spherical SNR 



S(v). 



£ 



P(v)r 2 dr . 



(16) 



Note that here we also take into account the geometrical di- 
lution effect that takes place in spherical shocks (in contrast 
to planar shocks). The highest energy particles were injected 
when the shock was young and, therefore, had a much smaller 
surface area than at the observation time. Potentially different 
injection rates at earlier times are involved through parameter 
b. 

Spectrum d!6l > may be written in a dimensionless form for 
the Sedov SNR: 



r 



da r^(a)r(a) B(a,G Q )K(a) 



(17) 



Jo Wo-bB 



' exp 



fi(a, ©o, a>, d, e/ef) 



de 



where r g = dr/da and e — v/v c (E m . dx , B a ). Note that, contrary 
to the case of uniform plasma, the electrons with energy e now 
have critical energy s c = e 2 <x B Z?. 

Spectrum d!7t depends on a set of parameters, O , oj, d 
(i.e. on B ), £f (i.e. on the combination B e ff,s£max0 and b, while 
the "fit"-spectrum (|6jl is only a function of a. It is clear from 
Eq. illl that the shape of the spectrum emitted by the electron 
population immediately postshock (i.e. at a = 1) is exactly the 
same for any set of these parameters, because B, K, and S equal 
unity immediately after the shock. The shape of this spectrum 
coincides with the shape of the "fit"-spectrum l|6} with a = 1 . 
The downstream evolution of the postshock spectrum depends 
on a set of parameters. 

Examples of the spectra at a — 0.8 are shown in Fig. [2] The 
function S is almost independent of d and obliquity, if consid- 
ered close to the shock (up to a * 0.8). Therefore, spectra for 



different d and O are quite close. In agreement with Reynolds 
( 1998|, the maximum losses are on maximum electron ener- 
gies. Thus, the end of the spectrum becomes steeper with time 
(with the decrease in a). Lines 1, 5, 6, 7 on Fig.[2]show spectra 
for different laws of evolution of K. It is apparent that no value 
of b produces the broadening of the spectrum. The only pos- 
sibility of reducing losses at high energies is the high value of 
6jf oc (B e ff,s£max0 1 (li ne 4). Lines 1, 2, and 3 demonstrate the 
role of ISM nonuniformity. The shape of the spectrum is almost 
insensitive to a> while the amplitude is sensitive. If the ISM 
density increases (co < 0), then the amplitude of the spectrum 
is less. If the density decreases, then plasma produces more 
emission at the same a. The stronger gradient of the post-shock 
density for smaller u is the reason for this behaviour. 



The spectra calculated with Eq. fl!7i are shown in Fig. [3] 
The higher the electron energy, the greater the radiative losses. 
Therefore, the superposition of spectra from different regions 
inside the volume of SNR might a priori be expected to be nar- 
rower than the initial postshock spectrum. For example, line 1 
may be approached by the fit-spectrum with a = 1.1; i.e. it is 
steeper than the initial post-shock spectrum at a = 1 . All the 
spectra presented in Fig.[3]lie between fit-spectra with a = 0.8 
and a - 1.3. Note that Fig.|3]shows spectra for those sets of in- 
put parameters that produce maximum differences in spectra at 
a = 0.8 (Fig. [5). None of these sets produces a broad spectrum 
(I17> . which could be approximated by Sf\ t with a « 0.5. The 
single case of the perpendicular shock (line 9) produces a bit 
broader emission spectrum (close to with a * 0.8) in com- 
parison with the parallel-shock case. The reason for such be- 
haviour is that electrons with the same energy e are responsible 
for the synchrotron emission with higher energy (e c = e 2 cr B B) 
in the perpendicular shock (cr B = 4 for O = 7r/2 while cr B = 1 
for © = 0). 



The efficiency of the space parameter variation in the mod- 
ification of the shape of the high-energy end of the synchrotron 
spectrum is low. For example, line 4 differs from the line 1 
maximum 1.25 times. The reason for this is a rapid fall in the 
amplitude of spectra with the decrease in a caused by syn- 
chrotron losses and hydrodynamic properties. Therefore, most 
of the contribution to the integral illl is given by the regions 
quite close to the shock, with a « 1. 



2.2. Time evolution of E n 



The model in Sect. l2.1l assumes that the shock accelerates elec- 
trons to the same maximum energy E m . dx at any time. It could 
be expected that the shock - during its evolution - is able to 
accelerate particles to different maximum energies. Reynolds 

oc yi with q 



( 1998 ) considered models where £ mal oc V q with q — 1 if the 
maximum energy is limited by radiative losses, or q = 2 if £m ax 
is determined by the finite time of acceleration. 
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log(e) log(e) 

Fig. 2. Evolution of the synchrotron spectrum downstream of the shock. The thick solid line shows spectra for different sets of 
parameters at a = 1. Other lines correspond to spectra at a = 0.8. The basic set of parameters is © Q = 0, a> — 0, d — 1, 6f = 1, 
b — 0. The corresponding spectrum is shown by line 1. Other lines show spectra with only one parameter changed: 2 - a> = -5, 
3-w=l,4-3 = 10,5-d = 0.1,6-6 = 2,7-6 = 5,8-6 = -2,9-© o =ff/2. 



If the maximum energy varies in time as E„ 
the spectrum (fl*6l may be calculated as 

f 2" 

S(s)cc I da r%(a)r(a) B(a,® )K(a) 

J<3min 



V(f) 9 , then 



e 2 exp 



(18) 



£(a, o ,w, of, e/ef) 

v/VcC^maxCO.Bo), 



Here we redefine e = E/E mlix {t) and e 
where f is the time of observations. 

In media with to < 3, positive q broadens spectra for differ- 
ent a compared with the case q — 0. The larger q, the broader 
the specrum for the given a. However, this broadening is too 
small to have a prominent effect on the integral spectrum Jl 81 . 
The time evolution of E mlix is inefficient in producing broader 
artificial spectra. Namely, for the basic set of parameters (see 
caption to Fig.|2j, the spectrum Jl 81 with q = 1 (q = 2), nor- 
malized at e = 10~ 4 , differs from the spectrum with q = 
maximum 1.15 (1.35) times. 

The reason for this inefficiency lies again in the rapid de- 
crease in B and K downstream of the shock that causes small 
weight in the integral Jl 8i of the downstream layers with 
broader spectra in comparison with the layers a » 1 . 




Fig. 3. High-energy end of the synchrotron spectrum fl!7i . 
flmin = 0.9. Lines correspond to models listed on Fig. [2] All the 
presented spectra are normalized at e = 10~ 4 . The fit-spectrum 
(|6j coinsides with line 1, if a = 1, and is shown by line 10 for 
a = 0.5. 



2.3. Variation in parameter on the surface of the SNR 

In this section we consider neither the downstream evolution 
of the initial spectrum generated by a shock nor the dis- 
tribution of parameters downstream. The role of the parameter 
variation inside the volume of SNR in the modification of the 
artificial spectrum is minor, as shown above. The purpose of 
the present section is to study a possible variation in the pa- 
rameters on the shell of SNR. Three of them can influence the 
synchrotron spectrum: B, K, E m . dx . Their changes on the shock 



can be caused by non uniformity of the ambient magnetic field 
(B) and/or obliquity (K, E mia ). In this paper we are consider- 
ing the only uniform magnetic field and spherical SNRs, hence, 
the nonuniformity of ISM where the density changes along the 
surface of the SNR is not studied here (it could be one more 
reason for the variation in K (which depends on density) and 
£ mal (which depends on shock velocity)). 

Let 9 be an angle between the direction of B Q and the given 
direction. The variations in parameters are given by E max (8) = 



6 
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Fig. 4. High-energy end of the synchrotron spectrum ( II 91 . 8 m - 
7r/6, a E = 1- Lines: 1 - 6k = °°, 0e — °°; 2 - 6k — oo, E — n/3; 
3 - 6 K = oo, 6 E = tt/6; 4 - 6 K = oo, 6 E = tt/9; 5 - K = tt/4, 
E = tt/9;6-6 k = 7r/9,0 E = tt/9;1-0 k = 7r/12,6fe = tt/9. AU 
presented spectra are normalized at e = 10~ 4 . The fit-spectrum 
© with a = 0.5 is shown by line 8. 

E t f E (0), K s (0) cc / K (fi), B s (fi) oc o-b(6»), where £» is a constant. 
Integration of the emissivity Q over the surface of SNR be- 
tween = -r 6 m (a polar cap) instead of dl7> gives: 



S(s) 



Jo 
Jo 



c/6- 



g- B (g)/ B (g)/K(e)tane 

cos 2 -yl - cos 2 m tan 2 



(19) 



eV B (6) 



e 2 exp 



de, 



where e = E/E m , e = v/v c (E* 
f E (0) = exp(a E (0/8 E ) 2) J 

/ K (c?)=e X p(-(0/0 K ) 2 ), 



B Q ). We adopt 



(20) 



(21) 



where E and 6k are some constants and a E is either 1 or -1. 
Expression \2\\ is able to approximately restore the depen- 
dence of the injection efficiency on obliquity from Ellison et 
al. i 



(1995 ) with6 K =tt/9h-^/4. 
Equation ( 1191 shows that the shape of the high-energy end 
of the emission spectrum does not depend on ,/k(6) if /e(0) = 
const (in this case the spectrum shape coincides with the shape 
of fit-spectrum S&t with a = 1). The gradients in the distri- 
bution of E max plays the main role in the modification of the 
shape of spectrum fl!9l >. If a E — 1, then one has a broader spec- 
trum (Fig.@] lines 1-4), and a narrower one if a E — -1. If E max 
varies, then changes in K also affect on the shape of the spec- 
trum (Fig.0] lines 4-7), but its role is less efficient. 



of the high-energy end of the synchrotron and inverse Compton 
spectra of SN 1006 and G347.3-0.5. This broadening consists 
in a < 1 and is an intrinsic property of processes involved in 
the formation of the spectrum in the region of acceleration, as 
in a model of Ellison et al. (2000 1. 

2. The broadening of the observed spectrum could be an 
artifact of observations. This means that electrons were accel- 
erated on the shock and their real a — 1, but the observed syn- 
chrotron spectrum looks as it is produced by electrons with 
a < 1. This could be caused by the fact that the observed 
spectrum is a superposition of spectra from different regions of 
SNR. However, both the variation in parameters downstream of 
the shock (even strong variation if SNR is in a nouniform ISM) 
and the time evolution of the injection efficiency and/or E m . dx 
are ineffective in the production of a broader artificial spec- 
trum. The shape of this a spectrum is in general narrower than 
the shape of the spectrum emitted by the electron population 
immediately after the shock. Such an artificial spectrum can 
be approximated by S at with ax 1 1.3 for different sets of 
parameters. Only perpendicular shocks are able to put the ob- 
served synchrotron spectrum close to S&t with a = 0.8 due to a 
higher compression of the magnetic field that makes the emit- 
ted photons harder (the critical frequency s c oc <x B (@ )). 

3. If the broadening of the spectrum is due to the variation 
of parameters, there has to be a variation in E max on the sur- 
face of SNR. The gradients dE m . dx /d0 have to be fairly strong 
to produce spectra close to Sfa with a = 0.5. Namely, if the 
spectrum is observed from a region of SNR, like a polar cap 
between = -r m , then E =s 20 m /3 and, therefore, E max has 
to change inside this cap almost 10 times. 
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3. Conclusions 



1 . The internal property of the electron spectrum produced on 
the shock can be one of the reasons for the observed broadening 



